On the possible secondary component of the order parameter observed in London penetration 

depth measurements 

A. Valli 

Dipartimento di Fisica, Universita Sapienza, P.le Aldo Mom 2, 1-00185, Roma, Italy and 

Institute fiir Festkorperphysik, Technische Universitat Wien, Wiedner Hauptstrasse 8//138, 1040 Vienna, Austria 

G. Sangiovanni 

Institute fiir Festkorperphysik, Technische Universitat Wien, Wiedner Hauptstrasse 8//138, 1040 Vienna, Austria 



o 



(N 

"a 
o 
o 

Oh: 



O 

o 



> 

in 

rn 

O 

o 



X 



M. Capone and C. Di Castro 
ISC-CNR and Dipartimento di Fisica, Universita Sapienza, P.le Aldo Moro 2, 1-00185, Roma, Italy 

We discuss the effect of a secondary component of tlie superconducting order parameter on tiie superfluid 
density in the cuprates. If we assume a main d^2 _y-2 gap, the most stable realization of a mixed order parameter 
has a time-reversal breaking d^2_y2 + idxy symmetry. In this state the nodes are removed and the temperature 
dependence of the superfluid density changes from the linear behavior of a pure d-wave to a more rounded shape 
at low temperature. The latter is compatible with the behavior experimentally observed in the in-plane magnetic 
field penetration depth of optimally doped La2-xSrxCu02 and YBa2Cui,OT-s- 

PACS numbers: 74.20.Rp, 74.20.Fg, 74.25.Dw 



I. Introduction 



The identification of the pairing mechanism behind high- 
temperature superconductivity in copper oxides ' remains one 
of the greatest challenges in solid state physics. A key in- 
gredient is the symmetry of the order parameter, which is ex- 
pected to reflect that of the pairing interaction thus provid- 
ing information on the microscopic mechanism. The well- 
established evidence of lines with vanishing amplitude in the 
gap function of cuprates along the V-X direction of the Bril- 
louin Zone (nodes) indicates a dominant cf^2_j,2 symmetry 
of the order parameter^, hardly compatible with the stan- 
dard phonon pairing mechanism, which leads to an isotropic 
s-wave order parameter. Anyway a small secondary com- 
ponent of the order parameter can develop either sponta- 
neously or driven by external factors like magnetic field, 
doping or presence of magnetic impurities^"^. The devel- 
opment of a mixed order parameter has been also invoked 
to explain anomalies observed in the thermal conductiv- 
ity in magnetic field of Bi2Sr2CaCu208^. Moreover sub- 
stantial deviation from the d^i^-y^-wave symmetry has been 
clearly observed in YBa2Cu307_5 (YBCO) both in tunnel- 
ing measurements^ and in laser angle-resolverd photoemis- 
sion spectra revealing nodeless bulk superconductivity*. A 
series of low-temperatures anomalies has been observed in the 
in-plane magnetic field penetration length in muon-spin rota- 
tion (^SR) experiments^ '^. Experiments in optimally doped 
La2_i:Sr2,Cu02 (LSCO) and YBCO have indeed shown a 
low-temperature bump superimposed to the linear tempera- 
ture behavior associated to d-wave superconductivity and to 
the presence of nodes. These deviations from d-wave behav- 
ior have been associated to a secondary component, which 
has been proposed to be isotropic s-wave in light of its vul- 
nerability to a magnetic field. Among alternative proposals, 
some'^''^ do not assume the presence of a mixed order pa- 
rameter, and they associate the low-temperature feature to a 



non-local response of the d-wave superconductor, which mod- 
ifies the magnetic field distribution in the vortex state with re- 
spect to the standard London model. Beside the interpretation 
in terms of a secondary superconducting s-wave'^, a particle- 
hole secondary gap associated to spin density wave ordering 
has been invoked'^. 

Here we focus on the secondary superconducting gap in- 
terpretation, and we show that only a d^2_y2 + idxy mixed 
order parameter can reasonably describe the jjlSR experimen- 
tal results. Previous analysis"* has shown that, assuming a 
leading dx2-y2 symmetry, the most stable realization of a 
mixed order parameter has indeed this time-reversal break- 
ing dx2-y2 + idxy symmetry. Moreover the development of 
such a time-reversal breaking order parameter does not require 
ad hoc assumptions, in contrast e.g. with a d + s symmetry, 
which requires completely unrealistic parameters as long as 



the dr 



is the dominant component of the order parameter. 



This work is organized as follows: In Sec. II we present 
our model and approach. In Sec. Ill we discuss the general 
behavior of the superfluid density with a mixed order param- 
eter and presents the comparison with experiments. Sec. IV 
contains our conclusions. 



n. Model 

In this section we briefly summarize the formalism used in 
Ref . 1 8 to identify the conditions for a secondary component 
to establish in the presence of a dominant dj;2^y2 wave. We 
consider a two-dimensional square lattice characterized by the 
€4^1, point group and a single band with dispersion 
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where t and t' are the nearest and next-nearest hopping param- 
eters, /i is the chemical potential and a = 1 is the lattice spac- 
ing. Values for hopping parameters for different compounds 



have been chosen according to density-functional theory cal- 
culations in the local-density approximation'^. 

The aim of the present analysis is the understanding of the 
competition between the different components of a supercon- 
ducting order parameter. Therefore we do not attempt a so- 
lution of a microscopic model including different kind of re- 
alistic interactions, and we simply consider an effective low- 
energy interaction, whose strength in each symmetry chan- 
nel controls the corresponding instability. Moreover, we will 
study the superconducting phase within the Bardeen-Copper- 
Schrieffer (BCS) mean-field approach, which fully takes into 
account for the symmetry of the order parameter This ap- 
proach is reasonably justified for instance by the relatively 
large doping of the samples of Refs. 9,12 

We now briefly recall some relevant aspects of the BCS 
equations for a mixed order parameter, referring to 18 and 
references therein for more details. If we require the in- 
variance under the symmetry of the lattice of the modulus 
of the order parameter, the latter has to transform either as 
an irreducible rappresentation or as a complex combination 
of the form A^ + ilSy (with A^ and A" transforming as 
two different irreducible representations) which breaks time- 
reversal invariance. The development of each harmonic with 
a given symmetry is controlled by a specific spatial compo- 
nent of the pair potential. The isotropic s-wave is associated 
to the local component of the potential Vq, which is repulsive 
in the cuprates due to the strong Coulomb interaction. The 
d^i_yi and extended-s {Sr^ij^yi) are controlled by the nearest- 
neighbor coupling Vx, while the d^y and Sj.y (which are anal- 
ogous to dx2-y2 and Sx2^y2 with lobes along the diagonal di- 
rections in the plane) are related to the next-neighbor coupling 
V2. Here we will simply assume that Vq is repulsive and that 
Vi and V2 are attractive. For the sake of definiteness we report 
the equations for the dj.2_y2 + id^y mixed order parameter 
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Here /3 = 1/T is the inverse temperature, ajd(k) ~ 
cosi\<ixa) — cos(kj^a) and ^^'(k) = 2sm(k^a)sm(k.ya) are 
the harmonics associated to dx2-y2 and dxy-v^sse respec- 
tively. Ad an d A^' are the associated com ponents of the gap 
and Ek = V^^ + A>^(k) + A^,cj^,(k). A^, A^- and the 
chemical potential are derived solving self-consistently Eqs. 
(2). An energy cutoff wq is used in the first two k-sums. 

For realistic dispersions, the dj.2_y2 symmetry is the lead- 
ing instability for small dopings due to the Van Hove singular- 
ity (VHS)^". When the main d.j.2_y2 order parameter appears 
at Tc, for T < Tc the effective dispersion ek is gapped and 
any secondary instability requires a minimum (critical) value 
for the associated interaction strength, as opposed to the case 
of an instability developing in a Fermi sea ground state. The 
dxy component turns out to be the best candidate for the sec- 



ondary gap (i.e., it has the lowest critical value of the inter- 
action) since it has the largest contributions from the regions 
in which the main gap has nodes. Since the onset of a sec- 
ondary component is essentially determined by the competi- 
tion with the main gap, one can favor a mixed state by re- 
ducing the dx2-y2 component. The complementarity between 
dxy and (i^2 _y2 also implies that the two components can exist 
symultaneously for a wide range of parameters. Other insta- 
bility channels require much larger couplings and, even more 
importantly, hardly give rise to a "coexistence" of order pa- 
rameters. In most cases, and in particular for s-wave com- 
ponents, the secondary order parameter can less efficiently 
exploit the Fermi-surface portions in which the first gap has 
nodes. Therefore, if we increase the associated coupling, we 
have an abrupt change from a pure dx2^y2 to a pure order 
parameter of different symmetry, and a very fine tuning is re- 
quired to have both order parameters. 

The focus of this paper is the effect of a secondary compo- 
nent of the superconducting order parameter on the superfluid 
density ps, which is directly related to the London penetra- 
tion depth by the relation A^^ = Aire^ps/mc^, being m the 
electron mass and c the speed of light, ps is defined as 
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where the first term is the zero-temperature contribution, 
while the other is the current-current response. For BCS pair- 
ing, in case of spin degeneracy, the previous expression then 
reads 
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(4) 
being ./(ck) = l/(e'^''' + 1) the Fermi distribution function 
for the Bogoljiubov quasiparticles. 



III. Superfluid Density in tlie mixed state 



Before addressing the comparison with experimental data, 
we consider the effect of the onset of the ^3.2 _ j,2 + id^y mixed 
order parameter in general terms. In Fig. la we plot the tem- 
perature dependence of the superfluid density and of the com- 
ponents of the superconducting gap (all normalized to their 
T = values) for Vi/t = 0.50, ¥2/1 = 1.05, loq/I = 0.25 
and t' /t ~ 0.25'"*. Indeed the results show that the linear 
temperature behavior characteristic of the d-wave state'^, as- 
sociated to the presence of nodal quasi particles, is modified 
below a temperature T'^ (see Fig. la). The low-temperature 
feature of the superfluid density is clearly related to the de- 
velopment of a dxy gap, which fills the nodes of the dx2-y2 
component below the secondary "critical temperature" T^. In 
this regime, in which a dx2^y2 + idxy order parameter is sta- 
ble, the shape of the superfluid density is more similar to that 
of a usual s— wave superconductor, reflecting the absence of 
low energy excitations. This shows that an "s-wave-like" be- 
havior at low temperatures does not automatically suggest an 



s-wave component, and that the dx2_y2 + id^y order param- 
eter generates a temperature behavior which reproduces the 
qualitative results of Refs. 9,10. 

We now briefly discuss how the shape of the superfluid den- 
sity depends the parameters of the system. A crucial parame- 
ter which varies in the different cuprates is the next-neighbor 
hopping t"^ which controls the position of the VHS. There- 
fore t' can push the singularity close to the chemical potential, 
thereby favoring the d^^-y-i at the expenses of the secondary 
gap. Indeed, as shown in Fig. lb) at fixed doping T'^ decreases 
as the chemical potential approaches the VHS. The same be- 
havior holds for the amplitude of the secondary gap as ex- 
pected within BCS. Similar results are naturally obtained by 
changing the hole concentration instead of t' , i.e. shifting the 
chemical potential and preventing it to lie within the cutoff 
energy range from the VHS. In practice the variation of t' in 
different materials can be quite large, and it affects the sym- 



metry of the order parameter much more than the doping, if 
the latter is taken in the physically relevant regime. In some 
cases (e.g. in LSCO compounds, where t'/i ~ 0.15''), the 
chemical potential can approach or cross the VHS in the rel- 
evant doping range (see Fig. Ic). On the other hand when t' 
is larger (e.g. in YBCO compounds) and the singularity is far 
from the Fermi level, the effect of doping becomes less impor- 
tant. Also orthorombic distortions or bilayer splitting reduce 
the dx2-.y2 gap, allowing for a larger secondary component'^. 
As a more technical note, the value of the cutoff cjq plays a 
role in the stability of the secondary component because it se- 
lects the portion of density of states which contributes to the 
effective coupling, i.e. a small cutoff makes the system more 
sensitive to the details of the bandstructure. For the range of 
parameters of interest this reflects in a stronger effect of the 
VHS, which favors the main component at the expenses of 
the secondary one. 
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FIG. 1: (a) ps and superconducting gaps A^2_y2 and Axy as a function of T. Each quantity is normalized to its T = value in order to 
confront the curves: it is clear as the low-temperature feature of the superfluid density relies on the existence of a secondary component of the 
order parameter (b) Behavior of ps for different values of t' /t for 5 — 0.17; (c) Behavior of ps for different dopings at fixed t' /t — 0.15 
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We now turn to the experimental evidences discussed above 
considering the specific cases of optimally doped LSCO' and 
YBCO'^. We use parameters (Vi/t = 0.55, Fa/t == 1.1, 
Wo A = 0.25, t'/t = 0.135 for LSCO and Vi/t ^ 1.1, V2/t == 
1.275, ujo/t = 0.25, t'/t = 0.35 for YBCO) that reproduce 
the experimental dispersions and the zero-temperature value 
of the gaps. The doping is 5 = 0.17 in both cases. As shown 
in Fig. 2, our simple theoretical approach well reproduces 
the temperature behavior of ps for a wide range of temper- 
ature. The appearance of the secondary component is much 
more pronounced for LSCO, in agreement with the above 
analysis about the role of t'/t. The deviation between the 
BCS results and the experiments close to Tc are obviously ex- 
pected because of the relevance of fluctuations for quasi two- 
dimensional strong-coupling superconductors. 

Our analysis shows that the experimental evidence of a low- 
temperature "bump" on top of the linear temperature depen- 
dence can be understood in terms of a dr^2_y2 + id^y order 
parameter without invoking an s-wave component. On the 
other hand an isotropic component is extremely hard to estab- 
lish when the main order parameter has d.j.2_y2 symmetry'^. 



In particular an s-wave order parameter requires a huge local 
attractive interaction, in clear contradiction with the almost 
universally recognized role of Coulomb repulsion. We em- 
phasize that even if we accept an attractive local component, 
the s-wave and the (ij.2_j,2-wave gaps hardly coexist in the 
same Fermi surface. It is therefore almost impossible to repro- 
duce the experimental behavior with an isotropic secondary 
gap without making totally ad hoc and unrealistic assump- 
tions and a very fine tuning of parameters . 

As a final remark we focus our attention on the effect of an 
external magnetic field, which seems to flatten out the low- 
temperature behavior of the penetration depth in the experi- 
mental data'"'". Several conflicting interpretations have been 
proposed. Some of them'-^'^ do not rely on the presence of a 
mixed order parameter, but relate the flattening to a non-local 
response of the d-wave superconductor, which modifies the 
magnetic field distribution in the vortex lattice with respect to 
the standard London model. Other studies identify the low- 
temperature feature with a second gap being either spin den- 
sity wave'^ or a different superconducting gap in the same 
spirit of the present analysis'^. In Ref 16 in particular, the 
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FIG. 2: Theoretical results for superfluid density against ex- 
perimental data for the (normalized) spin depolarization rate on 
-I'll. 855*7-0.15(71104' and YBa2Cu30e.95^^ compounds. 



hardly shed lights on the symmetry of the secondary compo- 
nent. Within the second gap interpretation, we notice that the 
main element for the stability of the secondary component is 
actually the size of the main gap, and that tiny variations of 
the latter may completely suppress the former. 



IV. Conclusions 



In this paper we have shown the effect of a time-reversal 
breaking order parameter dj.2_y2 + id^-y on the temperature 
evolution of the superfluid density within a BCS formalism. 
As previously shown, this combination is the most stable 
mixed order parameter if the main component has 0^2,2 _y2 
symmetry. Moreover, it is essentially the only way to have 
a smooth evolution from a pure d-wave to a superconduct- 
ing phase which displays a secondary component at low tem- 
perature. The same smooth evolution is mirrored in the tem- 
perature behavior of the superfluid density, in which a small 
"bump" is superimposed to the linear behavior characteristic 
of a pure (ij.2_j,2-wave. 

We compared numerical results to experimental data on 
two cuprates and showed that the low-temperature feature ob- 
served in /^iSR measurements can be reproduced assuming 
reasonable parameters for the system in such an unconven- 
tional symmetry phase. 



fragility of the secondary component to an external magnetic 
field has been advocated as a proof of its 's-wave nature. It is 
crucial to observe that the experimental magnetic fields (0.02, 
0.1, 0.64 T)*^ are too low to directly affect the secondary gap, 
whatever the mechanism could be. In this sense, the effect 
of the field can only be a minor indirect consequence, and it 
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